Interaction of two level systems in amorphous materials with arbitrary phonon fields 
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To describe the interaction of the two level systems (TLSs) of an amorphous solid with arbitrary 
strain fields, we introduce a generalization of the standard interaction Hamiltonian. In this new 
model, the interaction strength depends on the orientation of the TLS with respect to the strain 
field through a 6 x 6 symmetric tensor of deformation potential parameters, [R] . Taking into account 
the isotropy of the amorphous solid, we deduce that [R] has only two independent parameters. We 
show how these two parameters can be calculated from experimental data and we prove that for any 
amorphous bulk material the average coupling of TLSs with longitudinal phonons is always stronger 
than the average coupling with transversal phonons (in standard notations, 7; > 7t). 

PACS numbers: 63.50.+X, 61.43.Fs, 61.43. Er 



The thermal properties of dielectric crystals at low 
temperatures are well described by the Debye model. If 
the temperature is much smaller than the Debye tem- 
perature of the crystal, then the optical phonon modes 
are not excited and the only contribution to the heat 
capacity and heat conductivity comes from the acous- 
tic phonons. In three dimensional (3D) systems, low- 
frequency acoustic phonons have a linear dispersion rela- 
tion, lo = Ct,ik, where u> is the angular frequency, Cj and 
Ci are the transversal and longitudinal sound velocities, 
respectively, and k is the absolute value of the phonon's 
wavevector, k. This gives a specific heat proportional to 
the temperature to the power three (cy ocT 3 ). 

A good estimate for the heat conductivity, k, is k = 
|cycZ, where c is an average sound velocity and I is 
the phonon mean free path, which depends not only on 
the material, but also on the sample quality. Impurities 
or lattice defects, even at low concentration, reduce the 
phonon mean free path and in this way may decrease 
dramatically the heat conductivity— £££ As a result, the 
temperature dependence of the heat conductance is de- 
termined by the dependence of the phonon mean free 
path on its energy. Since such dependences can be very 
much different for different phonon scattering mecha- 
nisms, the resulting temperature dependence of k can 
be in general rather complicated. 

In high-quality crystals of relatively small size and at 
sufficiently low temperature, the phonon mean free path 
may become comparable to or bigger than the crystal di- 
mensions. In this case the phonons will scatter mainly at 
the surfaces and the mean free path is limited by the 
surface diffusivity and the geometrical features of the 
sampled In this case I is independent of the phonon fre- 
quency and k (x cy cx T 3 . 

Continuing to decrease the temperature and the size 



of the system, we get into the mesoscopic regime, where 
one or more dimensions of the system become compa- 
rable to the dominant phonons wavelength. Typically 
we find this regime in nanometer-size objects, at tem- 
peratures of a few Kelvins or less. At such scales, the 
phonon interaction with the surfaces becomes important, 
since it leads to coupling between different vibrational 
modes. This can lead to a very complicated set of phonon 
modes, with nonlinear dispersion relations. One exam- 
ple are the Lamb modest, which are among the eigen- 
modes of a free standing infinite membrane. Due to the 
boom of nanotechnology, mesoscopic systems are used 
in more and more practical applications (see for exam- 
ple Refs. 0@iIIlIllHmiH- The physical properties 
of mesoscopic devices can differ dramatically from those 
of bulk systems, the difference being more pronounced at 
lower temperatures. 

Important parts of many mesoscopic devices - such 
as microbolometers, electromechanical sensors and actu- 
ators - are ultrathin free standing membranes. The ther- 
modynamics and thermal transport of such membranes 
are specifically important at low temperatures, where 
heat release is a bottleneck for device performance. At 
the same time, the low-temperature transport properties 
are rather unusual. For example in many experiments the 
heat conductivity, k, of large, thin membranes or long, 
narrow and thin bridges is proportional to T p , where p 
takes values between 1.5 and 24 i 9 i 10 i 11 The specific heat 
of mesoscopic membranes is more difficult to measure. 
It can, however, be extracted from the amplitude of the 
temperature oscillations in AC measurements and it ap- 
pears also to be proportional to T p , where p lies roughly 
between 1 and 2J£ 

The temperature dependence of k and cy in meso- 
scopic membranes and bridges can partly be explained 



by a crossover from a three-dimensional to a two- 
dimensional phonon gas distribution. This crossover 
takes place when the dominant thermal phonon wave- 
length, Xt = iirhc/k^T, is comparable to the membrane 
thickness. 

Because of the finite thickness, there are gaps in the 
phonon spectra, the phonon dispersion relations near the 
band edges being nonlinear. As a result, the temperature 
dependence of the specific heat deviates from the law 
cy oc T 2 , which one would intuitively expect from a two- 
dimensional phonon gas. Instead, one obtains for very 
low temperatures cy oc T . 10 ^ 1 ' 13 ^ 5 

To study the heat conductivity in mesoscopic insu- 
lators, bridges have been cut out of the above men- 
tioned membranes. There, a decrease of the exponent 
of the temperature dependence of n from 2 to about 
1.5 was observed. Furthermore, the cut-off frequency of 
the temperature oscillations in AC heating measurements 
(f c oc n/cy) was observed to increase with temperature 
for the narrowest bridges . 10 ! 16 ' 17 Taking into account the 
fact that the edges of the bridges are very rough due 
to the cutting process, the measured data could be ex- 
plained by using the dispersion relations of the above 
mentioned Lamb modes to calculate heat capacity and 
heat conductivity^ 

Nevertheless, in some experiments the same behav- 
ior, k oc T p with p being roughly 2 ) 9 ' 10 i 14 ' 16 ' 17 and the 
increase of f c with the temperatur e 10 ' 16 was observed 
also above the 2D-3D crossover temperature. To explain 
these features we have to extend our model and take into 
account the amorphous structure of the material - low 
stress amorphous silicon nitride - used in the above men- 
tioned experiments. 



I. TWO LEVEL MODEL 

Amorphous or glassy materials differ significantly from 
crystals, especially in the low temperature range, where, 
for 3D bulks, n oc T 2 and cy oc T4J£ These temperature 
dependences were explained by the presence of specific 
dynamic defects. These defects are modeled by an en- 
semble of so-called two level systems (TLS) that exists in 
the material . 19 ' 20 A TLS can be understood as an atom, 
or a group of atoms, which can tunnel between two close 
minima in configuration space, forming a hybridized dou- 
blet state. The presence of such minima is a hallmark of 
the glassy state. If the energy splitting between these 
minima is < ^bT, then the TLS can be excited from its 
ground state onto the upper level. In this way the TLSs 
contribute to the heat capacity. TLSs can also scatter 
phonons and in this way decrease their mean free path 
and, correspondingly, the heat conductance. 

An effective double-well potential and the tunneling 
of the atom between the two wells are depicted in Fig. 
[I ] 19 ' 20 Written in the 2D Hilbert space spanned by the 
ground states of the two wells, the effective Hamiltonian 




FIG. 1: An atom of a group of atoms moves in an effective 
potential like this. The separation between the ground states 
in the two wells, A, is much smaller than the energy scale of 
the oscillation frequency in the wells, wo- 
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with A describing the tunneling between the two wells. 
In general, A is called the asymmetry of the potential and 
A is called the tunnel splitting. The Hamiltonian {1} may 
be diagonalizcd by an orthogonal transformation O, 



h tls = T H Ths O - -g z = 2(0 -e 



(2) 



where e = VA 2 + A 2 is the excitation energy of this TLS 
and by the superscript T we denote in general the trans- 
pose of a matrix. Since the TLS can be in only two states, 
let us denote the ground state by | J.) and the excited state 
by ||). Phonons at frequencies close to the level splitting 
e are strongly scattered by TLSs. 

The Hamiltonian parameters A and A are distributed 
with the density VP(A, A), where V is the volume of 
the solid. According to the standard tunneling model 
(STM), P(A, A) is assumed to have the form 



P(A,A) = P /A, 



(3) 



where Po is a constant. If expressed through the variables 
e and u = A/e, the distribution function is 



P(e,u) 



Pi 



(4) 



The strain caused by a phonon or any other kind of 
distortion of the material adds a perturbation to PtlSi 
which we denote by Hi. The total Hamiltonian, H = 
Ptls + Hi , enables us to describe the coupling between 
a TLS and the phonon field. In the STM the variation of 
the off-diagonal elements of Ptls are neglected, so Hi is 
diagonal^^* 2 ^ 
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The perturbation S is linear in the strain field, Si 3l — 
and in general may be written as 6 = 2jijSij. Here, 
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as everywhere in this paper, we assume summation over 
repeated indices. The 3x3 symmetric strain tensor is 
defined as £y = -^(diUj + djiii), with Ui, i = 1,2,3, be- 
ing the components of the displacement field. According 
to Eq. ([5]), a static strain would just slightly renormal- 
ize the eigenvalues of the Hamiltonian. Since thermal 
and transport properties of an insulating solid are deter- 
mined by lattice vibrations we are interested in the effect 
of a timc-dcpcndcnt strain entering the Hamiltonian Hi , 
describing the interaction of the TLS with a phonon field. 
In such a situation, the interaction produces transitions 
of the TLS between its eigenstates in the unperturbed 
state. These direct transitions contribute to the phonon 
mean free path and in this way influence the heat con- 
ductance. 

The effect of the time-dependent modulation of the 
spacing between the TLS energy levels is very important 
for low- frequency phonons 3i The populations of the TLS 
levels also change in time, but they lag behind the modu- 
lation of the interlcvel spacing. The resulting relaxation 
causes energy dissipation and, in its turn, phonon damp- 
ing. We do not consider this effect here since the heat 
conductance is governed by thermal phonons for which 
the resonant interaction is most important. 

In the STM the TLS interacts with 3D, transversally 
or longitudinally polarized, plane waves, which have all 
the components of the strain tensor proportional to the 
absolute value of the wave vector, k. The longitudinal 
wave will produce a compressional strain (S yz = S zx = 



S X z = S yx = 0) and the transversal wave 



will produce a shear strain (S xx = S yy = S zz = 0). 
In this case, the expression for 5 is always reduced to 
8 a = 2N a - r y a k (T , where a = t or I, denotes the transversal 
or longitudinal polarization of the wave, N a is the nor- 
malization constant of the phonon's displacement field 
and 7c is called the deformation potential parameter or 
coupling constant. Since N a has dimensions of length, 
N a k a is dimensionless and j a has dimensions of energy. 
To calculate the transition rates of the TLS from one 
cigenstate to another, induced by the interaction with a 
phonon, we have to write Hi in the basis that diagonal- 
izes H Ths : 



H[ = T HiO 
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The off-diagonal terms in H[ determine the transition 
rates. Note that, since the excitation and de-excitation of 
the TLS cause the absorption and emission of a phonon, 
respectively, 6 is implicitly of the form S oc S^b^ + h.c, 
where denotes a phonon annihilation operator. We 
omit this form here, but we shall use it explicitly in Sec- 
tion QIB] To calculate the transition rates, let us denote 
the population of phonon modes by n and the population 
of excited TLS states /. For example in thermal equilib- 
rium, nk CT = [e' 3r,Wk - - 1] _1 and f e = [e^ e + l] -1 , where 
(3 = l/(ksT). In these notations the TLS de-excitation 
amplitude into a phonon of wave-vector k and polariza- 



tion <7, due to phonon-TLS interaction, is 



(ma, + 1,1 \H[\n ka t) = -1 



hk 
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-7<r-vWr + 1- (7) 
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Using ([7]) , one obtains the contribution of a phonon with 
wave vector k and polarization a to the TLS transition 
probability due to phonon emission and absorption, re- 
spectively: 



r cm (e,k, cr) = 7, 
and 

r abs (e, k, cr) 



2 ^ A 2 , 
2 irk A 2 



+ l)S(hc a k - e) , (8) 



rikaSihc^k - e) . 



(9) 



In the relaxation time approximation, summing 
r a b s (e, k, a) — r em (e, k, cr) over all the phonon modes, we 
obtain the TLS relaxation time a o 18 ' 21 
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Similarly, we obtain the phonon relaxation time by sum- 
ming over all the TLS statesi 18 i 21 



irhujk 



hpc 2 a 



7 2 P • tanh 



\ 2 



(11) 



The low temperature mean free path of a phonon in 
the amorphous material, Cg-r^u, can be determined from 
the so-called unsaturated ultrasound attenuation, i. c., 
from the attenuation of an external acoustic wave of such 
small amplitude so that only a small fraction of TLSs are 
excited out of equilibrium. From the acoustic attenuation 
one can directly extract the product 7 2 -Po- Another way 
to determine 7 2 -Po experimentally is by measuring the 
relative shift in the sound velocity i 18 i 23 
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where To is a reference temperature at which Ac^ = 0. 
By these methods J. Black calculated the product 7 2 Po- 
In fused silica, for example, its values are within the in- 
terval (1.4 -j- 4.6) x 10 7 J/m 3 for longitudinal phonons 
and (0.63 — 0.89) x 10 7 J/m 3 for transversal phonons, see 
Ref. |23j and references therein. 

The values obtained for 7 2 i"b by these two independent 
methods can be combined to calculate the heat conduc- 
tivity as 



«(T) 
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which can be measured by yet another experiment.— 

Although sometimes there arc pronounced differences 
between the values for 7 2 Po obtained in different exper- 
iments (we are not concerned here with classification of 
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results) , it seems that always 7; > j t ■ To the best of our 
knowledge, this aspect is not explained in the literature. 
Moreover, the simplified expression for the perturbation 
term <5 in the interaction Hamiltonian does not allow us 
to calculate the interaction of the TLS with an arbitrary 
strain field. Since, as mentioned above, the elastic modes 
in mesoscopic systems have rather complicated displace- 
ment and strain fields, we need a more general interaction 
Hamiltonian which should incorporate the microscopic 
symmetry of the material around the TLS and the ori- 
entation of the TLS with respect to the strain field. In 
the next section we will build such an interaction Hamil- 
tonian and, after we apply it to 3D bulk systems, we will 
show how we can extract information about its parame- 
ters from experiments. The relationship 7; > 74 > is 
a natural result of our model. We will provide the re- 
sults for the interaction between TLSs and Lamb modes 
in ultra-thin membranes elsewhere. 



II. THE GENERAL TLS-PHONON 
INTERACTION HAMILTONIAN 

To generalize the TLS Hamiltonian we shall use the 
full expression for 6 of Eq. ([5]), 

S = 2^ ij S ij = 2[^:[S], (14) 

where 7^ are the components of the 3x3 tensor [7] and 
":" is the symbol for the dyadic product. Let us now find 
the general properties of [7]. 

As follows from Eq. (fl4j) . only the symmetric part of 
the tensor [7] has a physical meaning. Indeed, the dyadic 
product between a symmetric and an antisymmetric ten- 
sor is zero, so the antisymmetric part of [7], even if ex- 
istent, would not influence the results. We therefore as- 
sume [7] to be symmetric. Since all the tensors we have 
in our model are symmetric, it is more convenient to work 
using the abbreviated subscript notation, as described for 
example in Ref. 0. We will explain this method briefly. 

Let us assume that [A] is a symmetric 3x3 tensor (i.e. 
Aij = Aji). From its 9 elements, only 6 are independent. 
To get rid of the redundant 3 elements and also to make 
the tensors manipulation easier, we can write [A] in the 
form of a six component vector, A, in the following 2 
ways: 

A X = A XX , A 2 EE Ayy , A 3 ee A zz , 

A4 = A yz , Ac, = A zx , Aq = A xy (15) 

or 

Al EE A xx , A 2 EE Ayy, A 3 EE A zz , 

Ai ee 2A yz , A 5 ee 2A ZX , A 6 = 2A xy (16) 

Then the dyadic product of two symmetric tensors, [A] : 
[B] , may be written as [A] : [B] — A T ■ B if one of the 
tensors is written in the form (|15|) and the other one in 
the form (|16|) . Applying this way of writing to equation 



(O, we define 7 ee (>y xx , j yy , j zz , j yz ,j zx , j xy ) T and S ee 
{S xx , Syy, S zz /2S yz , 2S ZX , 2S xy , ), so that 5 = 27 • S. 

As in Ref. @, the subscripts in the abbreviated sub- 
script notations will be denoted by capital letters, 
/, J, K, . . ., which run from 1 to 6. 

A. The properties of the deformation potential 
tensor 

The vector 7 should characterize the TLS and its de- 
formability in the presence of a strain field. As explained 
before, the TLS is imagined as a particle or a group of 
particles that tunnels from one potential well to another. 
This tunneling may happen as a translation between the 
wells, or as a rotation^ In either case, there is a direc- 
tion associated to the TLS, which we call t, for example 
the direction defined by the two potential wells or the 
axis of rotation. One can expect that the orientation of 
the TLS (i. e., t) relative to the phonon's propagation 
direction and polarization has an effect on the interac- 
tion strength. The three components of t are the only 
co(ntra)variant quantities that describe the TLS, from a 
very general point of view, i.e. without building a micro- 
scopic model of the TLS. With these quantities we can 
construct the simplest symmetric 3x3 tensor 

(t X t X ty t X t Z \ 
t X ty ty tyt Z I t ' t , 

tx^z iytz ^ x J 

and a general one, [7] = [R] : [T] (i.e. 7 fc; = R^uT^), 
with Rijki = Rijik for any k and I. We shall also choose 
Rijki = Rjiki-, since the summation RijkiTij allows us 
to use this simplification. In abbreviated subscript nota- 
tions, T = {t 2 x ,t 2 y,tl,2tyt z ,2t z t x ,2t x ty) T , and R ijk i be- 
comes Rij in an obvious way. Then we can write 7 as 

7=[Rf-T. (17) 

Since T characterizes the orientation of the TLS, the rel- 
evant deformation potential parameters are contained in 
[R]. To make an analogy, the tensor [R] is similar to the 
tensor of elastic stiffness constants from elasticity theory. 
Still, the matrix [R] cannot be taken arbitrary. Like the 
clastic stiffness constants, the matrix [R] is determined 
by the local symmetry properties of the atomic lattice, 
around the TLS. We will deduce here the properties of 
[R]. 

We start by noting that the product hi = T T • [R] ■ S is 
a scalar, so it should be invariant under any rotation of 
coordinates. Moreover, since [R] reflects the microscopic 
symmetry of the lattice around the TLS, we shall do a 
sequence of transformations to obtain the properties of 
[R]. For the moment, let us assume that the lattice is 
simple cubic: 

• Let us choose Si = S 2 = S 3 = S 5 = S e = 0, S 4 ^ 
and t — x. Then we rotate the coordinate sys- 
tem through 7r about the z axis. Under this rota- 
tion [R] is invariant, whereas T = (i^,0,0,0,0,0) T 



5 



transforms into T' = T and S = (0, 0, 0, S 4 , 0, 0) T 
transforms into S' = -S. Then, h x = T • [R] ■ S = 
RlAt 2 x S A = T' ■ [R] ■ S' = -RutlSi implies that 
i?i4 = 0. Performing similar rotations about the 
other axes of the coordinate system with appropri- 
ately chosen T and S, we can show that Ru = 
for any I = 1, 2, 3 and J = 4, 5, 6. 

Taking S = (Si ± 0, 0, 0, 0, 0, 0) T and T = 
(0, t\, t 2 z ,2t y t z , 0, 0) T (so t x = and t v , t z ^ 0) and 
rotating the coordinate system through tt about the 
z axis, we transform S into S' = S and T into 
T' = (0,tl,t*,-2t y t z ,0,Q) T . From the product 



arrive at 



h[ — Si(R2%ty + Rzit z 
Rzit z + Ra2t y t z ) 



RAi2t v t z 



Sl(i?2lt 2 



hi, we deduce that Ru = 0. 



Again, performing similar rotations about the other 
coordinate axes we can show that Rjj = for any 
7=1,2,3 and J = 4,5,6. 

Up to now we proved that the matrix [R] is block- 
diagonal. Let us see if we can simplify it further. 

• Assume that t x — 0, whereas t y and t z are dif- 
ferent from zero. If, moreover, from the compo- 
nents of the strain vector only S5 is different from 
zero and we rotate the coordinates through n about 
the x axis, then T = (0, t 2 y , t 2 z , 2t y t z , 0, 0) T -> T = 
(0,t 2 y ,t 2 z ,2t y t z ,0,0) T and S = (0, 0, 0, 0, S 5 , 0) -> 
S' = (0,0,0,0,-5*5,0) Calculating the products 
T T .[R]-S and (T') T • [R] • S' wc find R 45 = 0. 

By similar arguments we conclude that -R45 = 
i?46 = R&6 = R54 = R§± = R&5 = 0. 

• Due to the cubic symmetry, h\ and H\ should not 
depend on the notation of axes. Therefore Rn = 

R22 = -R33) R12 = R13 = R23, R21 = R31 = R32, 

and i?44 = i?55 = i?66- 

These are all the constraints that we can impose on [R] if 
the lattice around the TLS has cubic symmetry. Now let 
us make one more simplification and assume the system is 
isotropic and find a relationship between the parameters 
R11, R12, R21, and -R44. 

Again from elasticity theory we know that both T 
and S transform under a rotation of coordinates like 
T' = [N] ■ T and S' = [N] ■ S, where the matrix [TV] 
is defined for example at page 75 of Ref. @. Since at a 
rotation of coordinates only T and S change, and not 
the deformation potential tensor (the interaction Hamil- 
tonian should look the same in any coordinate system), 
hi = T T ■ [R] ■ S = T T • [N} T ■ [R] ■ [N] ■ S for any T and S. 
This implies [R] = [N} T ■ [R] ■ [N]. Taking arbitrary ro- 
tations about each of the coordinate axes, we obtain the 
final conditions: R12 = R21 = R13 = R31 = R32 = R23 
and i?n — 2i?44 = R\2- 

By denoting Rn = 7, Ru/j = Cj an( i -R44/7 = ^, we 
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(18) 



with C + 2£ = 1. 

The form of [R] is very general, yet, the parameters 
7, £ and implicitly £ may vary from one type of TLS to 
another. For example, for the two types of TLSs men- 
tioned in the beginning of this subsection-translational 
and rotational (see Ref. [25j and references therein)- the 
parameters of [R] may be different. 



B. Physical results 

We can now assemble back the expression for 5: 

r]-S. (19) 



S = 2 7 T T • 1 
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a = (1 


S 



For calculations of physical quantities, wc have to write 
Hi in second quantization^ First, we introduce the ex- 
citation and de-excitation operators for the TLS, 



(20) 



These matrix operators obey Fermionic commutation re- 
lations and satisfy the conditions: a z = (2a^ — 1) and 
a x = (cv + a). The phonon creation and annihilation 
operators are denoted by and 6 M , respectively. Here 
we use n to denote general phonon modes. For 3D plane 
waves, /1 = (ker). 

With these definitions, in the basis that diagonalizcs 
Htls, the interaction Hamiltonian reads: 

7 A r 



Hi = ——T T . [r] • [SA + S*&t] (2a f a - 1) 

~^T T ■ [r] ■ J2 [SA + S*6+] (^ + a) , (21) 



where we used the notation Hi to distinguish this form 
from the one in Eq. ([6]). By S M we denote the strain 
produced by the phonon field at the position of the 
TLS. 

In first order perturbation theory, the matrix element 
for the absorption of a phonon (kcr) by an unexcited TLS 
is 



(n kCT ,T \Hi\n ka + 1,1} = 



jA / hrik a 



2Vpu) ka 



[r] ■ S 



kcr ■ 



(22) 

Applying Fermi's golden rule, we obtain the transition 
probability for this process. Assuming that in an amor- 
phous solid the TLS directions are uniformly distributed, 
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wc sum the contributions of the TLSs averaging over their 
directions and we obtain an average transition rate, 



r i"k CT ,T)> kCT + l,i) ~ C <r7 



Vpc a € 



—5(huj ka -e), (23) 



where C a is a constant that depends on the polarization 
of the phonon. Similarly, the emission rate is 

^ ~2 ( n k<7 + l)7rfc A 2 
1 |n kCT +l,T),l«k„,i) = L al y~. ^J d (^k CT - ej • 

(24) 

Equations (f2"3T) and (f2"4T) are similar to equations © and 
([5]) respectively, with 7 2 replaced by Ctj 2 and 7 ; 2 replaced 
by C7 2 . The constants Ct and C are 

C; = -L(15-40£ + 32£ 2 ) (25a) 
15 

Ct = (25b) 

Even though we are not able to make any statement 
about the range in which £ takes value, we can still make 
the interesting prediction: 



Ci> C t > for any £, 



(26) 



in agreement with the experimental data. 

The TLS relaxation time and phonon absorption time 
can still be calculated by Eqs. (fT0|) and (fTT|) . with 7 2 
and 7 2 replaced by C7 2 and C7 2 . Therefore, C7 2 and 
C7 2 can be calculated from unsaturated ultrasonic at- 
tenuation or sound velocity shift experiments. Once the 
values C7 2 and C/7 2 are obtained, the ratio 



4C 
C t 



15 

e 



40 

T 



32 



(27) 



gives us the value of £, which further enables us to cal- 
culate (. If the value of 7 can be extracted from phonon 
echo experiments^ then all the elements of the displace- 
ment potential tensor [R] are known. 

As a numerical example, we take from Ref. [23| the two 
sets of values for Po7; 2 an d P)7 2 that we used before. 
For the first one^ where Po7 ; 2 = 1-4 x 10~ 5 J/m 3 and 
P 7 2 = 0.63 x 10~ 5 J/m 3 , the two solutions for £ are 
£1 = 0.55 and £ 2 = 1-2. Using the formula P07 2 = 
P 7 2 /C t (£) (see Eq. [25b| we obtain (P 7 2 )i = 7.8 x 10 7 
and (Po7 2 )2 = 1-7 x 10 7 . For the second set of values 2 ^- 
P07 2 = 2.0 x 10~ 5 J/m 3 and P 7 f 2 = 0.89 x 10~ 5 J/m 3 - 
we get £1 = 0.55 and £2 = 1-2, which correspond to 



(Po7 2 ) 



11 x 10 7 and (P 7 2 ) 2 = 2.4 x 10 7 . Note 



that although P07 2 (or P07 2 ) changes significantly from 
one experiment to another, the ratio Po7 ; 2 /Po7 2 does not 
change much, which leads to close values for the param- 
eter £. From these measurements only, we cannot make 
the difference between £1 and £ 2 - 



III. CONCLUSIONS 

To describe the interaction of a TLS with an arbitrary 
strain field, we introduced a generalization of the stan- 
dard TLS-plane wave interaction Hamiltonian. Such a 
generalization is useful in the description of mesoscopic 
systems, where the phonon modes are not anymore the 
simple, transversally and longitudinally polarized plane 
waves, but have more complicated displacement fields 
and dispersion relations. We used the symmetry prop- 
erties of the system to deduce the properties of our in- 
teraction Hamiltonian, Hi. 

We showed that if the TLS is in an isotropic medium, 
then H\ depends on four free parameters: two from the 
orientation of the TLS, t, and another two, denoted 7 
and £ (fTB|) . which describe the coupling of the oriented 
TLS with the strain field. Since, in an amorphous solid, 
t has random, uniformly distributed directions, the first 
two parameters are averaged and the effective form of the 
Hamiltonian has only two free parameters, 7 and £, like 
the Hamiltonian of the standard tunneling model. From 
sound absorption, sound velocity change, heat conductiv- 
ity and phonon-echo experiments one can calculate these 
two parameters, like we did in the end of the previous 
section. 

The other way around, if for example the distribution 
of the TLS orientations is not uniform, knowing the val- 
ues of the parameters 7 and £ one can eventually recon- 
struct the distribution. 

Different types of TLSs-translational and rotational 
ones-may have different deformation potential parame- 
ters. For one type of TLSs we proved that 7; > jt> 
which is confirmed by all the experimental results we 
know about. 
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